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Abstract. We report H- and A"' -band speckle interferometric observations of the HAeBe star V376 Cas. Our observations show 
that the object is partially resolved. The visibility curves suggest three separate components: a large scattering envelope visible 
only in the H band, a component approximately 100 mas in radius, and a component with a Gaussian HWHM of approximately 
8±3 mas, corresponding to approximately 5 AU at 600 pc distance, which contributes most of the flux. We compare the smallest 
structure to the radius of dust sublimation in the radiation field of the star and find that the radius is approximately six times 
larger than that expected. This may indicate that the inner regions of the system are in fact obscured by a flaring circumstellar 
disk or torus seen close to edge-on. 
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1. Introduction 

Many intermediate-mass young stellar objects (YSOs) are sur- 
rounded by large quantities of circumstellar material, giving 
rise to infrared excesses which often dominate their spectra. 
Much of this material is often distributed in a circumstellar 
disk, through which accretion onto the central object may con- 
tinue. Whilst simple disk models originally explained many of 
the features of YSOs, there remained some outstanding prob- 
lems with a thin disk picture. Principally, for both flat pas- 
sive reprocessing disks and active accretion disks, AF^ oc A^"^^^ 
(Lynden-Bell & Pringle 1974; Adams et al. 1987) whereas the 
observed spectral energy distributions (SEDs) are in general 
flatter than this. This implies that either a modification of the 
disk model or the presence of some other component to the 
circumstellar material is necessary. 

Hillenbrand et al. (1992) examined the SEDs of a large 
sample of HAeBe stars and divided them into three groups; 
those which could be modelled with a geometrically thin disk, 
sometimes requiring an optically thin inner hole, those with 
a rising SED at long wavelengths, which could be modelled 
only with an envelope, and those which lacked either disks or 
envelopes. These categories were suggested as an evolution- 
ary scheme for circumstellar environments of HAeBe stars. 
Hartmann et al. (1993) suggested instead that a combination 
of disks and envelopes was generally necessary. Since then, 
disk-Henvelope models have been further investigated by sev- 
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eral authors, including for example Natta et al. (1993), Pezzuto 
et al. (1997), and Miroshnichenko et al. (1997, 1999). 

The presence of disks together with envelopes is strongly 
suggested by the observations of Mannings & Sargent (1997, 
2000) who studied several intermediate mass YSOs in the mil- 
limetre continuum and various molecular lines. For several 
sources, they found line profiles consistent with Keplerian ro- 
tation. In some cases the rotation curve could be spatially re- 
solved. The quantity of material inferred from millimetre con- 
tinuum fluxes was generally not consistent with measured ex- 
tinctions and a spherically symmetric distribution. 

The main alternative to disk+envelope models is some type 
of flared disk (e.g. Chiang & Goldreich, 1997). For a reprocess- 
ing disk, such models modify the temperature power law and 
flatten the SED at far-infrared (FIR) wavelengths. Malbet et al. 
(2001) and Lachaume et al. (2003) considered flared-disk mod- 
els in which the disk is heated by both stellar light and viscous 
dissipation and had success in fitting SEDs. However, simple 
flared-disk models cannot reproduce the near-infrared (NIR) 
'bump' between about 1 and 3 jum observed in many HAeBe 
systems. Duflemond et al. (2001) and Natta et al. (2001) there- 
fore proposed a new variant of the flared-disk model, intro- 
ducing an inner hole and thereby producing a bright inner rim 
which contributes more NIR flux than a simple flared disk. This 
inner rim can puff up and partially or completely shadow the 
more distant disk. Beyond any shadowed region, the outermost 
regions of the disk are the same as envisaged in the simple 
flared-disk picture, maintaining the flux at FIR wavelengths. 

In recent years, advances in speckle observing techniques, 
the advent of adaptive optics, and NIR long-baseline interfer- 
ometry have begun to allow insights into the spatial distribu- 
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tion of the circumstellar material and outflows associated with 
YSOs. For example, Leinert et al. (2001) used a speckle tech- 
nique to observe 3 1 HAeBe stars and found many to have en- 
velopes. The complex environment of the HAeBe star R Mon 
was studied by Close et al. (1997) with adaptive optics and 
by Weigelt et al. (2002a) with bispectrum speckle interferom- 
etry, revealing several structures related to the disk and out- 
flow. Similar bispectrum speckle observations of the high-mass 
YSOs S140 IRS 1 by Schertl et al. (2000) and Weigelt et 
al. (2002b) and of S140 IRS 3 by Preibisch et al. (2001) re- 
vealed complex outflows associated with these objects. Monin 
& Bouvier (2000) used an adaptive optics system with a resolu- 
tion of about 70 mas in the J band to image an edge-on circum- 
stellar disk in the young triple system HV Tau. Danchi et al. 
(2001) imaged the the massive young star MWC 349 A using 
aperture-masking interferometry and found what appeared to 
be an edge-on disk. This object was also observed by Hofmann 
et al. (2002) who used bispectrum speckle interferometry and 
also saw the same probable edge-on disk. Other edge-on disks 
have been imaged using HST (see for example McCaughrean 
& O'Dell 1996 or Brandner et al. 2000). Tuthifl et al. (2002) 
used an aperture masking interferometry technique to resolve 
the massive HAeBe star LkHo' 101 and saw a circular ob- 
ject about 50 mas across, corresponding to 3.4 AU at 160 pc. 
This object showed an asymmetry reminiscent of an almost 
face-on but marginally inclined disk. Malbet et al. (1998) re- 
solved FU Ori with the Palomar Testbed Interferometer (PTI). 
Akeson et al. (2000) observed four YSOs with the PTI and re- 
solved three of them, two T Tauri stars and an HAeBe star. 
Millan-Gabet et al. (2001) made interferometric observations 
of a sample of fifteen HAeBe stars and saw many resolved 
sources. Typical sizes were a few milliarcseconds, correspond- 
ing to between about 0.5 and a few AU. Some of the above 
results were surprising, since the hot regions of classical accre- 
tion disk models, from where the NIR flux should largely arise, 
should be too small to be resolved at this level. 

V376 Cas is a HAeBe star embedded in the small dark 
cloud L1265 and associated with considerable nebulosity. At 
optical wavelengths, on arcsecond scales, this nebulosity is 
seen to be elongated along an axis with position angle approx- 
imately 120° (Leinert et al. 1991; Piirola et al. 1992; Corcoran 
et al. 1995; Asselin et al. 1996). The source displays an ex- 
treme degree of linear polarization in the red, ~ 20 - 25%, 
higher than for any other intermediate mass YSO, indicating 
that a high proportion of the emerging light has been scattered. 
Polarization maps were produced at 900 nm by Leinert et al. 
(1991), in the / band by Piirola et al. (1992), and in the optical 
by Asselin et al. (1996). These maps show strong centrosym- 
metric polarization in the lobes of the elongated nebulosity and 
areas of low polarization situated near the star on either side of 
the extended envelope. Leinert et al. (1991) used speckle po- 
larimetry at K to show that this pattern persists to subarcsecond 
scales. This pattern of polarization suggests a bipolar reflec- 
tion nebula aligned northwest-southeast, which in turn implies 
the presence of an edge-on disk at a position angle of approx- 
imately 30°. HST WFPC2 images obtained from the archive 
(Fig D confirm this elongated source. However, at this higher 
resolution, the overall nebulosity is not reminiscent of an edge- 



on disk with position angle 30°, but rather suggests an out- 
flow cavity to the west of the object. This in turn suggests a 
nearly edge-on disk at position angle approximately 140°, to- 
gether with the nearside lobe of a bipolar outflow extending to 
the west. The observed light probably arises mostly from scat- 
tering from the inner edge of the outflow cavity or from the 
surface of a flared disk or circumstellar torus. 

2. Observations and data reduction 

The speckle interferograms were recorded in September 2002 
with the SAO 6 m telescope in Russia. The detector of the 
speckle camera was a Rockwell HAWAII array. The two fil- 
ters used corresponded approximately to H (central wavelength 
1.648 yum, bandwidth 0.317 /urn) and K' (central wavelength 
2.1 15 yum, bandwidth 0.214 //m). Further observational param- 
eters are listed in Table 1 . 

The modulus of the object Fourier transform (visibility) 
was obtained with the speckle interferometry method (Labeyrie 
1970). Speckle interferograms of calibrator stars, which were 
measured to be unresolved, were recorded just before and just 
after the object and served as reference stars for the determina- 
tion of the speckle transfer function (see TableQ. 

3. Results and discussion 

The reconstructed visibilities of V376 Cas are presented in 
Fig. |2] The upper two panels show the two-dimensional vis- 
ibilities in K' (left) and H (right). The visibilities show some 
evidence of elongation along a position angle approximately 
70° from the vertical. The //-band visibility is dominated by a 
central spike, indicating the extended emission in the 3.8" field 
of view of H, which is not evident at K' . The azimuthally av- 
eraged visibilities are shown in the two lower panels of Fig.|2l 
These plots have both been truncated at twelve cycles per arc- 
second, beyond which the quality of the data is significantly 
degraded by noise. There is clear evidence of resolution at both 
wavelengths. The /T'-band visibilities decline by about 20% 
and the //-band visibilities by about 30% over the range of spa- 
tial frequencies sampled. Both curves show a significant de- 
crease in steepness at ~4 cycles/arcsec. This structure suggests 
a two-component visibility model for the K' band and a three- 
component visibility model for the //band. The best-fitting 3- 
and 2-component Gaussian models for H and K' are shown in 
the plots. A second model was also fitted, in which the most 
compact Gaussian component (on the sky) is replaced by a 
uniform-brightness ring. The inner radius of the ring was the 
fitted parameter, and the outer radius was set at 110% of this. 
The intention of fitting this model is to allow direct comparison 
with the results reported for similar objects in the literature (e.g. 
Millan-Gabet et al. 2001). To determine uncertainties in the fit- 
ted parameters, each individual target-reference star pair was 
fitted separately and the resulting parameters averaged. There 
were five such independent observations in total in each band. 
The peak-to-peak variation of each parameter was then used as 
an estimate of the uncertainty. Table |2llists the parameters and 
associated uncertainties of all these fits. 
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Fig.l. HST WFPC2 image of 
V376 Cas taken with the F814W 
fiher at 800 nm. This image, which 
was obtained from the archive, is the 
combination of one 100 s and one 
300 s exposure, both of which were 
made in March 2000. North is up and 
East is to the left. The brightest knot 
probably corresponds to the component 
observed at H and K' . 



Table 1. Observational parameters. Ac and A/l give the central wavelength and bandwidth of the filter FoV is the field of view 
of both the data recording and data processing. A^r and A^r are the total numbers of speckle interferograms of V376 Cas and the 
reference star, respectively, and A^p is the number of target-object pairs observed per night. T is the exposure time per frame and 
S is the Seeing (FWHM). In the last column, the names of the reference stars are given. 
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The Gaussian fit components are an approximately 1" 
HWHM extended nebula visible only in H, an approximately 
100 mas HWHM component, corresponding to about 70 AU in 
radius at an assumed distance of 600 pc (Chavarria-K 1985), 
and a small, approximately 8+3 mas HWHM component com- 
prising a large fraction of the flux in each band, corresponding 
to a region of about 4.8 AU in radius at 600 pc (see Fig.|2j. 
It should be stressed that smaller components could in fact be 
present which are not resolved at our longest baselines. We can- 
not, of course, exclude the possibility of another break in the 
visibihty curve at baselines longer than 6 m. 

The H-band visibility function shows resolved structure 
both along and across the elongation. This, together with the 
approximately circularly symmetric appearance of the K'-band 
visibihty function allows us to rule out a binary model as the 
sole explanation for the resolved source and to identify at least 
part of the resolved structure with cuxumstellar material. It is 
possible that a companion is also present, with a position an- 
gle of approximately 160° or 340°, and contributes partially 
to the resolved structure. It is also possible that the elongation 



is due to the presence of a nearly edge-on disk. If this were 
so, the disk plane would run along a direction perpendicular to 
the observed elongation in the visibility function, that is, from 
northwest to southeast. Such a disk would be broadly compat- 
ible with the interpretation offered above for the HST image, 
showing an outflow cavity extending to the southwest. 

3.1. The 8 mas component: Dust sublimation? 

An explanation often found in the literature for the observed 
NIR sizes of YSOs is that the observed structure lies at the 
dust sublimation radius, that is, the radius at which circum- 
stellar dust sublimates due to the radiation field of the cen- 
tral star. This gives rise either to a hot inner edge of a disk, 
as in the model of Dullemond et al. and similar models, or to 
a sharp inner boundary to a spherically symmetric shell. This 
type of explanation was first suggested for individual objects 
by Tuthill et al. (2001) and Natta et al. (2001), and was inves- 
tigated for a sample of resolved YSOs by Monnier & Millan- 
Gabet (2002). Monnier & MiUan-Gabet (2002) demonstrated 
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Fig. 2. Upper panels: Two-dimensional visibilities of V376 Cas. Left column filter 2115/214nm. Right column filter 
1648/3 17 nm. The contours are drawn at levels 0.95, 0.90, 0.85, 0.80, 0.75 and 0.70. Lower panels: Azimuthally averaged visi- 
bilites. Left column filter 2115/214mn. Right column filter 1648/317imi. The individual dots show the observed visibiUties, the 
solid line the fitted model consisting of two Gaussians at 21 15 nm and three Gaussians at 1648 nm. 
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that the NIR sizes of various HAeBe and T Tauri stars appear- 
ing in the literature correlates with the stellar luminosity in a 
way that is consistent with their being determined by dust sub- 
limation. Here, we compare the smallest resolved component 
of V376 Cas with the expected dust sublimation radius for this 
object. 

The dust subUmation radius can be estimated as 

= l-l>^(io55^)'''(T^rAU, 

where Qr is the ratio of the incident absorption efficiency to 
the re-emitted absorption efficiency (Monnier & Millan-Gabet 
2002). 



For hot stars, Qr can vary strongly if the grain size is 
smaller than about one micron. For grains of 1 /im or larger, 
2r ~ 1, but Qr may reach values of around 50 or more for 
grain sizes of less than 0.25 fim and stellar surface tempera- 
tures of 10,000 K or greater However, if large grains existed in 
sufficient quantities, they would survive longer and determine 
the position of This logic was followed by both Monnier 
& Millan-Gabet (2002) and Tuthill et al. (2001), who both as- 
sumed Qr = 1. 



The maximum temperature at which dust grains can survive 
is estimated to be around 1750 K for graphite and about 1400 K 
for silicate grains (Hillenbrand et al. 1992). We use a sublima- 
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Table 2. Results of circular symmetric fits to the observed visibilities. Two different models were fitted at each wavelength. The 
first model consisted of either two (at K) or three (at H) Gaussians. For the second model, the broadest Gaussian in the visibilities 
(corresponding to the narrowest component on the sky) was replaced with a uniform-brightness ring, the inner radius of which 
was a fitted parameter and the outer radius of which was set to 1 10% of the inner radius. 
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Uniform-brightness ring plus Gaussian models 
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stellar Luminosity 

Fig. 3. The sizes of YSOs as a function of stellar luminos- 
ity, based on Fig. 1 of Monnier & Millan-Gabet (2002). The 
data points show the inner ring radii of uniform-ring fits to re- 
solved sources from Millan-Gabet et al. (2001), the measured 
size of LkHa 101 from Tuthill et al. (2002), and the HWHM 
for Gaussian fits to two classical T Tauri stars from Akeson et 
al. (2000). The cross marks the position of V376 Cas. The lu- 
minosity has been taken as 517 L© with an assumed error of 
a factor of two. The radius plotted is the inner ring radius of 
the smallest component fitted to the K' data. The straight lines 
mark the model with Qr - I for three dust sublimation tem- 
peratures, 1000 K, 1500 K, and 2000 K. The dashed line is a 
model with a grain size of 0. 1 fim and - 1500 K. 

tion temperature of 1500 K but also consider temperatures of 
1000 K and 2000 K. 



Various estimates of the luminosity of V376 Cas are avail- 
able in the literature. Chavarria-K. (1985) estimated the lumi- 
nosity of V376 Cas to be 517 Lq, based on his own extinction 
and distance estimate and the photometry of Allen (1972) and 
Cohen (1974). Abraham et al. (2000) observed V376 in the 
mid to far infrared with ISOPHOT and derived a luminosity 
between 1 and 1000 microns of 430 Lq. Since the vast ma- 
jority of the flux emerges in this wavelength range, this value 
is probably not a serious underestimate due to the wavelength 
truncation at either end. Natta et al. (1992) estimated the lu- 
minosity of V376 Cas to be about 150 Lq, based on a combi- 
nation of new FIR observations and existing optical and NIR 
measurements. For the purposes of this analysis, we adopt the 
luminosity of 517 given by Abraham et al. (2000), together 
with an uncertainty of a factor of two. 

With these assumptions (L, =517 Lq, T, = 1500 K, Qr = 
1) the dust sublimation radius for V376 Cas given by 
Equation H] is 0.79 AU. This rises to 1.78 AU if the dust is 
destroyed at a temperature of 1000 K. These values lie ap- 
proximately a factor of three to six below the HWHM of the 
smallest fitted Gaussian for V376 Cas, at about 4.8 AU. This 
is illustrated in Fig. |3l which is based on Fig. 1 of Monnier 
& Millan-Gabet (2002). This plot shows the relationship be- 
tween stellar luminosity and dust sublimation radius for dust 
grains of 1 jjm or larger and three diff'erent dust sublimation 
temperatures, as well as the dust sublimation radius for a dust 
destruction temperature of 1500 K and a grain size of 0.1 fim. 
The position of V376 Cas in the K' band is marked. The er- 
ror bars in the size are determined from the fitting procedure, 
as described in Sect.|3]above. The luminosity is assumed to be 
uncertain by a factor of two. No account is taken of possible 
errors in the distance of 600 pc. This is because an error in the 
distance would move the point parallel to the model lines in 
the figure and would therefore have no effect on any compar- 
ison between the data point and the models. Also shown are 



6 



K.W. Smith et al.: Speckle interferometry of the HAeBe star V376 Cas 



points for a number of other HAeBe and T Tauri stars taken 
from the Uterature (see caption for references). 



3.1 .1 . Comparison with models 

The most compact component of V376 Cas which we resolve 
is about six times larger than the basic dust sublimation model 
predicts and also appears large compared to the population of 
similar objects shown in Fig. 13 Two possible explanations for 
this are that the dust grain population is dominated by grains 
significantly smaller than 1 //m, and that the inner regions of the 
system are obscured by a disk or torus viewed close to edge-on. 

In Fig.|3lwe show the estimated dust sublimation radius for 
grains of 0.1 /im in size. Such a grain size leads to a predicted 
dust sublimation radius significantly larger than the measured 
size of V376 Cas, demonstrating that even a modest change in 
grain properties might accommodate our result. The scattering 
envelope seen at H but not at K' may indicate that the scatter- 
ing cross section is very much smaller at K' than at H. This 
could imply that the population of dust grains in this extended 
region is dominated by small grains. However, the position of 
V376 Cas in Fig.|3lis not only somewhat high compared to the 
models but also compared to other similar objects. Invoking a 
small dust-grain size to explain this would therefore involve 
suggesting that V376 Cas is in this respect an exceptional case. 

Can the luminosity estimate for V376 Cas be in error? The 
spectral type of V376 Cas was determined by Cohen & Kuhi 
(1979) from low resolution spectra to be B5. Earlier, Herbig 
(1960) estimated the spectral type of V376 Cas to be late B 
or early A. Rodriguez & Canto (1983) observed the field at 
6 cm and saw no corresponding radio source down to a level 
of approximately 1 mJy, which implies that the spectral type 
cannot be much earlier than mid B. Comparison with evolu- 
tionary tracks (e.g. Palla & Stahler 1993; Driebe, 2001) sug- 
gests that a B5 star of surface temperature 14,000 K should 
have luminosity of between approximately 200 and 1000 Lq. 
This is broadly consistent with the Chavarria-K. (1985) and 
Abraham et al. (2000) estimates. However, if the source has a 
disk, the escaping radiation field will be anisotropic, and if the 
disk is viewed close to edge-on, this would lead to the photo- 
metrically determined luminosity, and hence the expected dust 
sublimation radius, being underestimated, probably by a factor 
of approximately two to four. Such an underestimate is easily 
compatible with a small uncertainty in the spectral type. 

Finally, if the NIR flux arises mostly at the inner edge of 
a disk, as in the model of Lachaume et al. (2003), DuUemond 
et al. (2001), and Natta et al. (2001), this region might wefl 
be obscured from view by the outer parts of any flaring disk 
or circumstellar torus. The observed flux would then be re- 
flected from the circumstellar material, and the size of the dom- 
inant component in the NIR might well be considerably larger 
than the original dust sublimation radius. This possibility is of 
course entirely consistent with the high level of linear polariza- 
tion observed from the object. 



3.2. The other components 

The largest structure seen at //but not at K' clearly coiTesponds 
to the structure seen by Leinert et al. (1991) and also by Piirola 
et al. (1992) and others. This extended //-band component is 
caused by scattered light from the extended nebula. This in- 
termediate structure (200 mas ~ 120 AU or 60 AU in radius) 
may represent scattering from the flaring parts of the outer disk 
and probably coiTesponds to the bright knot seen in the HST 
images. 

4. Conclusions 

Our H- and /T'-band speckle interferometric observations of 
V376 Cas suggest two resolved compact components, with 
characteristic radii of 100 mas and 8 mas. Additionally, a scat- 
tering envelope about 2" in diameter is seen in the //band. The 
resolved components together account for at least 25% of the 
flux at K' and 35% of the flux at H. The smallest resolved struc- 
ture may be related to the radius of dust sublimation at the inner 
part of a disk or envelope heated directly by stellar radiation. 
However, the actual size of this region (8 mas or 5 AU), is about 
six times too large for the supposed luminosity of V376 Cas. 
If V376 Cas possesses a disk, it is likely to be viewed close 
to edge-on, and this would lead to the photometrically deter- 
mined luminosity being an underestimate of the true luminos- 
ity. Furthermore, the outer part of any disk probably obscures 
the inner disk edge, where the dust sublimation region would 
be situated. The NIR flux observed would then most likely arise 
from scattering from the inner edge of a cavity above the disk 
plane. This hypothesis is further strengthened by the high de- 
gree of linear polarization displayed by this object. 
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